Origin of Low Thermal Conductivity in Nuclear Fuels 
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Using a novel many-body approach, we report lattice dynamical properties of UO2 and PUO2 and 
uncover various contributions to their thermal conductivities. Via calculated Griineisen constants, 
we show that only longitudinal acoustic modes having large phonon group velocities are efficient heat 
carriers. Despite the fact that some optical modes also show their velocities which are extremely 
large, they do not participate in the heat transfer due to their unusual anharmonicity. Ways to 
improve thermal conductivity in these materials are discussed. 



Today's nuclear fuels are based on 235 U and 239 Pu ele- 
ments where in a typical set-up, a nuclear reaction heats 
up a pellet made of either UO2 or its mixture with PuC>2 
and the heat is transformed to electrical energy. One of 
the major issues is to conduct the heat from a core of 
the pellet to its outer area which makes the evaluation 
of a high temperature thermal conductivity a key prob- 
lem. Unfortunately, the thermal conductivity of UO2 is 
very low and as a result, even a gradual increase of its 
value may lead to significant breakthrough in the perfor- 
mance of commercial nuclear reactors currently operating 
worldwide. 

In the present work, we argue that it is the unexpect- 
edly large anharmonicity of optical modes that results 
in a very low thermal conductivity of modern nuclear 
fuels. Consider semiconducting UO2 which is a main el- 
ement of UOX fuel, or a blend of U and Pu oxides which 
is a major substance of MOX fuel. The heat from the 
core of the pellet in these insulating systems is carried 
by phonons which are known to be very inefficient heat 
conductors. This brings a whole set of complex problems 
such as causing fuel pellets to crack and degrade prema- 
turely, necessitating replacement before the fuel has been 
depleted. 

Unfortunately, studying thermal conductivity [l[ as 
well as such properties as crystal structures, phase dia- 
grams, lattice dynamics and structural instabilities of the 
actinide based materials is a formidable theoretical prob- 
lem. Most of the previous works have concentrated on 
molecular dynamic simulations with empirically adjusted 
interatomic potentials However, the 5/ electrons in 

actinides are close to a localization-delocalization or Mott 
transition as it has been recently demonstrated for Pluto- 
nium metal Hdd). U0 2 and Pu0 2 are Mott-Hubbard 
insulators with energy gaps at both low and high tem- 
peratures. Despite an impressive set of past theoretical 
studies 0] 01 M [HI , the spectral functions of the high 
temperature paramagnetic regime cannot be obtained by 
static-mean field theories such as the Density Functional 
Theory [l2| and require a genuine many-body treatment. 

In the present work we use a novel electronic structure 
method [l3( capable of describing Mott insulating materi- 
als in order to address the structural properties and ther- 
mal conductivity of the modern nuclear fuels. Both UO2 



and PuG<2 are calculated using a combination of local 
density approximation (LDA) 12[ and Dynamical Mean 
Field Theory (DMFT) Q [l5( where relativistic 5/ shells 
of Uranium and Plutonium atoms are treated by exact 
diagonalization of corresponding many-body Hamiltoni- 
ans obtained by allowing a hybridization between the 5/- 
electrons and the nearest oxygen 2p orbitals. 

It is well known that a strong spin-orbit coupling of 
about leV present in actinides splits 14-fold degenerate 
/ level onto / 5 / 2 and fj/ 2 states. Group theoretical con- 
siderations assume that under cubic crystal symmetry, 
the fr,/2 6- fold degenerate level is further split onto Tg 
quadruplet and doublet. In both UO2 and Pu02, the 
Tg level comes approximately Q.leV below the T7 state, 
and valence arguments make it occupied by two electrons 
for the case of UO2 and fully occupied by four electrons 
for the case of PuO-2. This sequence of the levels dictates 
the low temperature properties of these two materials: 
The UO2 becomes magnetic where the 9-fold degenerate 
many-body ground state of the atomic /-shell 3 H4 with 
J = 4 is split onto 4 subsets, among which the lowest 
triplet state of T§ symmetry carries the moment of about 
1.7 — 1.8/xb below the Neel temperature of 30.8K flrjj ]. 
On the other hand, Pu02 has the filled Tg one-electron 
level, making its atomic 3 K4 multiplet split in such way 
that the non-magnetic Ti singlet comes lowest. 

The cluster exact diagonalizations which are carried 
out in our LDA+DMFT calculations support this atomic 
physics picture and only marginally alter the self-energies 
for the /-electrons from their corresponding atomic val- 
ues. In fact, the crystal structures of both materials show 
that the actinide elements are centered in the cubic en- 
vironment (see Fig.l) assuming that only relativistic Ti 
orbital (its shape is shown in Fig.l) of the / 5 / 2 state 
is properly coordinated by and strongly hybridized with 
the O 2p states. The corresponding hybridization func- 
tions A Q (w) are well fit by the single pole approximation, 
A a (u>) — \V a \ /(uj — P a ), producing only significant val- 
ues of the matrix elements V a ~ leV for in the j = 5/2 
manifold and several times smaller values for the Tg or- 
bitals. Similar effect is seen for Tg, T-?, Tg states in the 
j = 7/2 manifold. The position of the dispersive oxygen 
band centered at the pole P a comes out — 1 -. — 5eV be- 
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FIG. 1: Crystal structure of UO2 and PuCh: Uranium or Plu- 
tonium atoms (their shapes are shown by the relativistic 
orbitals) are centered inside the cubic lattice made by oxy- 
gen atoms. The arrows show the displacements of oxygens 
associated with the dispersive longitudinal optical mode. 



low the occupied f states, making both materials classical 
Mott-Hubbard insulators. We find that the fundamental 
energy gaps produced by the lower (occupied) and upper 
(unoccupied) Hubbard bands are equal to 2.2eV for UO2 
and 2.5eV for PuC>2. The experimental energy gap in 
UO2 is of the order of 2eV [l6j]. These data are obtained 
utilizing a recently developed matrix expansion algorithm 
[ri| which helps to perform full self-consistency with re- 
spect to both the charge densities as required by the LDA 
procedure and the 5 /-electron spectral functions as re- 
quired by the DMFT. Here, the 5/-electron self-energies 
extracted from the cluster exact diagonalization are sub- 
sequently fit using three-pole interpolation. An effective 
parameter U e ff = 3eV describing the on-site Coulomb 
repulsion among the 5/ electrons is used while the other 
Slater integrals (F^ 2 \ and F^) are computed from 
atomic physics, and are subsequently rescaled to 80% of 
their values to account for the effect of screening [l8j ]. 

In order to calculate the phonon spectra of UO2 and 
PuC>2 we utilize a new relativistic linear response phonon 
method which has been recently generalized for the 
LDA+DMFT scheme [lj} and has recently demonstrated 
its predicted power by finding the phonon spectrum of 5- 
Pu metal [2(| prior to the experiment 21 1 . In this linear 
response calculation the pole interpolation for the self- 
energies reduces the calculation of the dynamical matrix 
to standard linear response theory. Since our interest 
lies in the thermal conductivity at operating tempera- 
tures of nuclear fuels, we perform the calculations of 
the phonon dispersions for paramagnetic phases of the 
materials with the f-electron self-energies extracted at 
T = 1000.ftr , where not only the many-body ground states 



FIG. 2: Calculated phonon dispersions (filled symbols con- 
nected by lines) and experimentally measured phonons (open 
symbols) HH for UQ 2 . 



(r5 for UO2 and Ti for PUO2) but also various low-lying 
excited states are beginning to contribute. 

Fig. 2 shows the comparison between our calculated 
(filled symbols connected by lines) and measured (open 
symbols) [l6| phonon spectrum for UO2 along two sym- 
metry lines of the Brillouin Zone. Fig. 3 presents our pre- 
dicted phonon spectrum for PuC>2. Overall, both mate- 
rials demonstrate very similar dispersions. Due to var- 
ious numerical inaccuracies and approximate treatment 
of many-body effects, which, in particular, involves the 
use of U e f f , the overall mismatch between our theory and 
experiment for UO2 is about 15%. The same accuracy 
can be assumed for our predicted phonon dispersions in 
Pu02- Increasing the U e /f to 4 — 5 eV marginally affects 
the overall phonon spectra although makes all modes 
somewhat harder. 

Our calculations reveal the following features of lat- 
tice dynamics in these materials. First, acoustic modes 
are located at the low frequency region of the spectrum 
producing the group velocities of about 700to/s for trans- 
verse acoustic (TA) and llOOm/s for longitudinal acous- 
tic (LA) modes (see also Table 1). There are two groups 
of transverse and longitudinal optical modes which we 
schematically label as TOl, LOl and T02, L02. The 
first group represents lattice vibrations associated with 
the out-of-phase displacements of oxygens and little in- 
volvement of U or Pu atoms. The second group rep- 
resents the lattice vibrations of oxygens which vibrate 
in-phase against U or Pu atoms. 

There is a remarkably dispersive mode LOl which has 
a very large group velocity along (001) direction in both 
materials. This mode is primarily oxygen based and at 
the BZ center it is described by polarization vectors as 
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FIG. 3: Calculated phonon dispersions (filled symbols con- 
nected by lines) for PuCh. No experiment available. 

shown on Fig.l. We see that our theory correctly pre- 
dicts a softening of this mode as it approaches the zone 
boundary X point although it underestimates the exper- 
imental value of the frequency by a factor of two for UO2 
and, most likely, also for PuC>2. Note that similar dis- 
crepancies have been detected by us earlier in studying 
lattice vibrations of Plutonium metal [2(| . 

As we have gained access to the phonon dispersions, we 
are ready to discuss the lattice thermal conductivity in 
these systems which can be expressed for a cubic solid via 
the phonon frequencies w qj - , group velocities Vqj , phonon 
mean free paths Iqj and the Bose-Einstein distribution 
functions Nqj as follows [l|: 

\ - , dNqj 

qj 

It is generally known that various scattering processes 
determine the temperature dependent phonon mean free 
paths. These, for example, include normal phonon- 
phonon interactions where one phonon with wave vec- 
tor q gets scattered with the creation of two phonons 
with wave vectors qi and q2 such that q = qi + q2 , and 
Umklapp phonon scattering processes where the quasi- 
momentum conservation involves reciprocal lattice vec- 
tor G such that q = qi + q 2 + G. At high temperatures 
which we are interested in thepresent work, only Umk- 
lapp processes are significant [1], and the phonon mean 
free paths can be evaluated from the knowledge of the 
third-order anharmonicity coefficients, which can be ap- 
proximated by the dimensionlcss Griineisen constants de- 
scribing the change in the phonon frequency with respect 
to the atomic volume, i.e. jqj = din Uqj/dhi V. Hence, 
the efficiency of each phonon mode in the heat conduc- 



tion is directly related to its group velocity but reversely 
proportional to the square of its Griineisen constant 

Looking at the phonon spectra presented in Fig. 2 and 
Fig. 3, it is now clear that the best heat carriers are either 
the acoustic branches or the longitudinal optical branch 
LOl which, in fact, has anomalously large group veloc- 
ity along (001) (experimentally, its v^j is twice the LA 
branch!). One may then pose a question that why UO2 
or Pu02 are known to be such inefficient thermal con- 
ductors. When compared to most semiconducting solids 
where optical modes exhibit rather weak wave vector dis- 
persions and thus do not participate in the heat transfer, 
here the situation is much more favorable. The answer 
lies in the anomalously large anharmonicity associated 
with the LOl mode, making its mean free path signifi- 
cantly shorter. 

To make a comparative analysis, we have calculated 
the Griineisen constants 7^ associated with each vibra- 
tional mode at the T and X points of the Brillouin Zone. 
The results of these calculations are summarized in Table 
1 . As one can see the transverse acoustic modes have rel- 
atively large anharmonicity characterized by 7 = —1.41 
and 7 = —1.59 for UO2 and PUO2 respectively. These 
values remarkably decrease for the LA modes ( 7 = —0.50 
and —0.54) which at the same time show large group ve- 
locity. However, they become huge for the most disper- 
sive LOl mode. Here, the frequency at T shows positive 
7 (equal to 0.43 and 0.27) meaning that it decreases upon 
compression while the frequency at X shows negative 7 
(equal to —2.17 and —2.49) meaning its increase upon 
compression. Effectively, these two effects will add up to 
each other and result in the effective Griineisen constants 
equal to -2.60 for U0 2 and -2.76 for Pu0 2 . 

With these data in hands we are able to estimate con- 
tributions to the lattice thermal conductivity k from var- 
ious phonon modes. Using a simple Debye-like linear ap- 
proximation for the phonon dispersion and including cor- 
rections found empirically for the optical branches [22j ]. 
the results of such estimates at a temperature T = 1000 A' 
are presented in the last column of Table 1. We can con- 
clude that the only efficient heat carriers in both nuclear 
materials are their longitudinal acoustic phonons. The 
TA modes have relatively small contribution because of 
their smaller group velocities and relatively large 7, while 
the LOl modes also do not contribute due to its huge 
Griineisen constants which completely compensate the 
effect of its largest Vqj. The other modes have a negligi- 
ble influence on k. 

The total value of k in our calculation is a factor of two 
lower than the experimentally known thermal conductiv- 
ity of U0 2 equal to SSWm-iR- 1 at T = lOOOif [H]. 
This is naturally connected to the approximate charac- 
ter of treating anharmonic effects. However, we believe 
that the relative contributions from various branches are 
captured correctly in our estimates, which gives a funda- 
mentally new insight into the processes and factors that 
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TABLE I: Calculated phonon frequencies (in THz) and 
Griineisen constants at T and X points of the Brillouin zone, 
the phonon group velocities (xl0 2 m/s) estimated along (001) 
direction as well as contributions to lattice thermal conduc- 
tivity (ffm-'JT 1 ) at T = W00K for U0 2 and Pu0 2 . 

IKh 



control thermal conductivity in these materials. As a 
result, we finally discuss the origins of the large anhar- 
monicity of the most intriguing longitudinal optical mode 
LOl. As shown on Fig.l, this mode is primarily oxygen 
driven. Oxygens are arranged in a cubic lattice with Ura- 
nium or Plutonium atoms occupying every other center 
of the cube. It is then clear that the structure is far from 
close packing and the displacements of oxygens would 
involve large third order anharmonic effects. To fix this 
problem, one may think of mixing these materials with 
elements filling in the cubic interstitials of the lattice and 
preventing large ionic excursions. One example could be 
oxygen overdoped materials U0 2+x and Pu0 2+x where 
the latter one has been recently studied theoretically 0] • 
To summarize, using the LDA+DMFT method and 
linear response theory, we have studied lattice dynami- 
cal properties of UO2 and Pu02- Contributions to lattice 
thermal conductivity from various phonon modes were 
uncovered using the calculated group velocities and the 
Griineisen constants. It was found that the dispersive 
longitudinal optical modes do not participate in the heat 
transfer due to their large anharmonicity. Material de- 
sign of systems with the last effect supressed would open 
new possibilities to build more efficient fuels for modern 



nuclear industries. 
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